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The recent high-energy QCD equations including the effects of fluctuations are solved nu- 
merically. We discuss their asymptotic properties and compare them with the corresponding 
behaviour without fluctuations as well as with the statistical-physics model obtained in the 
saddle point approximation. The potential consequences on phenomenology, particularly on 
' ^ geometric scaling violations, are drawn. 
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In our way to understand the high-energy limit of QCD, a large amount of effort is made to 
obtain the evolution equations including the saturation effects. Indeed, when the system becomes 
\ more and more dense, one has to consider multiple scattering and recombinations, eventually 

leading to saturation and unitarity. The simplest of these equations is the Balitsky-Kovchegov 
(BK) equation^ which is known to admit traveling- waves solutions leading to geometric scaling 
Oh, T(k 2 ,Y) = T(k 2 /Q 2 (Y)), where Q S (Y) is the energy-dependent saturation scale. 

C^h, Recently, it has been proven 1 ^ that one also has to deal with fluctuation effects which become 

important in the dilute regime. At large N c , the evolution of the amplitude takes the form of a 
hierarchy of equations where is coupled to through BFKL evolution, to T^ k+1 ^ through 
merging leading to saturation, and to T^ k ~ 1 ^ through fluctuation effects coming from splitting. 

Interestingly enough, after a coarse-graining in impact parameter, this infinite hierarchy 
turns out to be equivalent to a Langevin equation which is formally equivalent to the BK 
equation supplemented with a noise term reproducing gluon number fluctuations 



d Y T(k) = a X (-d log(k) )T(k) - aT 2 (k) + a^/2Ka 2 T{k)u{k 1 Y), (1) 
where a = a s N c /i, x(l) 1S the BFKL kernel, and the noise u(k, Y) satisfies 

(v(k,Y))=0 and (v{k,Y)v(k' ,Y')) = —5(Y - Y')k5(k - k'). 

an 

The unknown fudge factor n fixes the strength of the noise and is expected to be of order 1. From 
this equation, one can see in which limits the different contributions to the evolution become 
important. On the one hand, the saturation contribution becomes relevant when T ~ T 2 i.e. 
when T ~ 1 which is the dense regime as expected. On the other hand, the noise term is 
important when T ~ ^a 2 T(k) or, equivalently, T ~ a 2 . This corresponds to the dilute regime 
where fluctuations are indeed expected to play a crucial role. 

In what follows, we p resen t the results'"' of numerical simulations of the Langevin equation 
The basic idea (see^EI f or a detailed description) for the numerics is to introduce the 
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(a) 100 events and their average (b) Dispersion as a function of Y 

Figure 1: Results of the event-by-event study. Y goes from to 50. 

probability P(5Y;Tq — > T) that the effect of the noise changes the amplitude from To to T in a 
rapidity step 5Y. P satisfies the Focker-Planck equation 

_ 2 

dvP(Y- T -» T) = [TP(Y; T -> T)] , (2) 
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which can be solved with the initial condition P(0;To — > T) = 5{T — To) and the effect of 
the noise can be computed at each spatial point using the associated cumulative probability 
Ft 0i sy(T). We perform one step of the evolution using 

Tnoise(^) = ^T(k,Y),SY ' 

T(k, Y + 5Y) = T noisc (k, Y) + SY [a X (-d L )T noisc (k, Y) - aT n 2 oisc (A:, Y)) , 

where u(k) is a uniformly-distributed random variable between and 1. For convenience, we 
shall use a = 0.2, L = log(fc) and R = «• The initial condition T(L, Y = 0) for the evolution 
is exp(— 2L) when L > and 1 — L otherwise. 

Before going into the analysis, let us note that if one performs a diffusive approximation 
in the BFKL kernel, Eq. (^Q) simplifies to the well-known stochastic FKPP equation. This 
approximation is known to give relevant results for the tail of the wavefront and allows one to 
study analytically^ the limit na^ — * 0. 

To begin with let us analyse a set of events generated through our method. As shown in 
Fig. ^ one clearly sees that each event shows a traveling-wave pattern, like in the BK case. 
In addition, the different events are spread. This dispersion means that although geometric 
scaling is preserved for each singular event, it is violated for the averaged amplitude. A more 
careful analysis (see the right plot of Fig. ^) shows that the dispersion, computed over 10000 
events, behaves like y/Y which is typical of these random processes. From the analogy with the 
stochastic FKPP equation, well studied in statistical physics, we know that, for asymptotically 
small values of no£. the dispersion should behave like log -3 (/-ca 2 ). However, in order to observe 
that behaviour, one needs values of kc? s as small as ^40~ 20 which cannot be reached by our 
method (see^for a model), It is important to note that, in the early stages of the evolution, we 
do not observe much dispersion, which may explain the fact that geometric scaling is observed 
in the HERA data. 

Another property that can be studied is the asymptotic speed of the traveling wave which 
is related to the saturation scale Q 2 S {Y) = kQexp(vY). For asymptotically small na 2 s , should 



decrease w.r.t. BK according to 
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Although we have not studied this asymptotic regime, we confirm this slowing down b as shown 
on Fig. |21 for realistic values of the noise strength. Note also that one has to evolve up to 
unphysical values of Y to observe the asymptotic velocity. 

Let us now proceed with the study the 
properties of the average amplitude, obtained 
from a set of 10000 events. The effect of the 
noise w.r.t. BK evolution without fluctuation 
is manifest in Fig. |3 As rapidity increases, 
the scaling violations become more and more 
important and the average front has a smaller 
slope. One also remarks that, although the 
asymptotic speed is smaller than for BK, in 
the early stages of the evolution, this speed ap- 
pears to be larger. For values of Y ~ 40, one 
sees however from the dashed curve in Fig. 01 
that this speed is already below the BK speed. 

Finally, we have compared the results of our 
simulations with the predictions of the stochas- 
tic FKPP equation, as shown in Fig. |1J On 
the left plot, presenting the results for (T) at 
Y = 20, it is obvious that the diffusive approx- 
imation (stochastic FKPP model) agrees very 
well with the numerical simulations c . The scal- 
ing violations for 10~ 4 < T < 1 are well ob- 
served and the rapid fall-off seen in the far tail 
of the front is coming from the fact that each 
event has a compact front. Theoretically, it can 
be seen from the probability P and, in physical 
terms, it means that the effect of the fluctua- 
tions is to cut off the ultraviolet tail. 

From a set of events generated using Eq. we can also generated higher-order correlators. 
The diffusive approximation predicts 

(T) « <T 2 > « • • • « (T k ) « . . . 

for asymptotic values of Y. As shown on Fig. 21 (right plot), numerical simulations qualitatively 
agree with this result for the tail of the front although the asymptotic regime is not yet fully 
reached. In the saturated domain, the diffusive approximation is not sufficient since it misses 
the logarithmic behaviour in the infrared. However, this region is characterised by (T 2 ^ ~ (T) 2 
as confirmed by simulations. 

It is interesting to notice that the sFKPP model agrees already at values of Y which are not 
extremely high as compared to the time needed to reach the asymptotic speed. Actually, it has 

'The asymptotic formula predicts a negative speed for c? s k, ~ 1 which proves the importance of subleading 
corrections. 

c The discrepancy at T > 1CP 9 is a numerical artefact. 



i 

0.98 
0.96 
Sg 0.94 
= >■ 0.92 
0.9 
0.88 
0.86 



~ i 1 1 r 



10000 15000 
Y 



1 

0.1 
0.01 
0.001 
0(BK 



20000 25000 



Figure 2: Asymptotic speed of the traveling wave for 
different strengths of the noise, compared with the BK 
asymptotic speed. 
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Figure 3: Evolution of the average amplitude (T) with 
rapidity Y = 0, 10, 50. One can see the effect of the 
noise by comparison with the BK curve. 
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(a) Form of the wavefront (Y = 20) 



(b) Numerical results for (T 2 ). 



Figure 4: Comparison of the numerical simulations (h — 1) with the stochastic FKPP predictions. Left plot: 
form of the wavefront. Right plot: correlations (T 2 \ together with IT) 2 and (T) for comparison. 



been proven^, in the asymptotically-small kq: 2 limit, that the formation time for the wavefront 
is of order /3 _1 7~ 2 log 2 (Ka 2 ) with (3 ~ 48.2. This is faster than the equivalent time for the BK 
equation and might have consequences on geometric scaling. Although this may depend on our 
choice of initial condition, it seems that one has a competition between a short time needed 
to reach the asymptotic form of the wavefront and the scaling violations induced by dispersion 
which (see Fig. ^) are not observed in the early stages of the evolution. This may leave a region 
in rapidity in which geometric scaling is verified. 

To summarise, we have seen from numerical simulation of the Langevin equation (^Q) that 
the contribution from fluctuations have a large impact on the high-energy evolution in QCD. 
It modifies the rapidity-dependence of the saturation scale and introduces geometric scaling 
violations. Our results, obtained for physical values of the coupling, are in agreement with the 
predictions coming from the stochastic FKPP equation. 

In the future, further analysis should be performed in order to clarify the question of geomet- 
ric scaling in the experimentally accessible range of rapidity where preasymptotic effects cannot 
be neglected. Furthermore, an extension of the solution to the complete hierarchy of equations, 
including the impact-parameter dependence, would be a precious part of our understanding of 
QCD at high-energy. 
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